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ABSTRACT 

We present the results of two XMM-Newton observations of the ultraluminous X-ray source 
(ULX) NGC 5204 X-1. The EPIC spectra are well-fit by the standard spectral model of a 
black-hole X-ray binary, comprising a soft multi-colour disc blackbody component plus a 
harder power-law continuum. The cool {kTin ^ 0.2 keV) inner-disc temperature required by 
this model favours the presence of an intermediate-mass black hole (IMBH) in this system, 
though we highlight a possible anomaly in the slope of the power-law continuum in such fits. 
We discuss the interpretation of this and other, non-standard spectral modelling of the data. 
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1 INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are luminous (Lx> 
10'^^ erg s~^), extra-nuclear X-ray sources in nearby galaxies. 
Their nature has been the subject of much recent debate (see Miller 
& Colbert 2004 for a review). In brief, the extreme X-ray lumi- 
nosities of many ULXs cannot be explained by stellar-mass black 
holes radiating isotropically at (or close to) their Eddington limit, 
leading to the suggestion that ULXs harbour a new 10^ - 10" 
Mq "intermediate-mass" class of black holes (IMBHs; Colbert & 
Mushotzky 1999). This hypothesis has gained strength from the re- 
cent detection of the soft X-ray emission characteristic of "cool" 
accretion discs, the signature of accretion onto an IMBH, in the X- 
ray spectra of several ULXs (e.g. Miller et al. 2003; Miller, Fabian 
& Miller 2004). However, there is strong circumstantial evidence - 
in particular the association of large populations of ULXs with very 
active star formation - that suggests that many ULXs may instead 
be a type of stellar-mass X-ray binary that somehow exceeds its Ed- 
dington limit (King 2004). Models for achieving super-Eddington 
luminosities from these systems include the "slim" disc model (e.g. 
Ebisawa et el. 2003), radiation-presure dominated discs (Begelman 
2002), mildly-anisotropic radiation patterns from sources emitting 
at or below the Eddington limit (King et al. 2001) and relativistic 
beaming (Kording, Falcke &. Markoff 2002). 

The best testing grounds for our models of ULXs are the near- 
est and brightest examples of this class of object. One such source 
is NGC 5204 X-1. It has a typical X-ray luminosity of the order 

2 - 6 X 10^' erg (0.5 - 8 keV), and is located ~ 15 arcsec 
from the center of a nearby (d — 4.8 Mpc) Magellanic-type galaxy 
(Roberts & Warwick 2000). It is most notable for possessing a stel- 
lar optical counterpart identified from its accurate Chandra posi- 
tion (Roberts et al. 2001), which HST observations resolved into 
two sources plus a third candidate counterpart (Goad et al. 2002). 



Further HST UV spectroscopy has now identified HST-3 as the 
ULX counterpart. Its spectral type is similar to a BO lb supergiant, 
but with peculiarities consistent with some Galactic X-ray binaries 
(Liu, Bregman & Seitzer 2004). Chandra showed its X-ray spec- 
trum to be a steep power-law continuum (with F ~ 2.8 — 3), and 
confirmed that it has been persistently X-ray luminous for 20 years, 
with flux variations of no more than a factor ~ 3 (Roberts et al. 
2004). The BO lb counterpart confirms that this ULX is a high- 
mass X-ray binary, though there are currently no constraints upon 
the mass of the putative black hole in this system. We present new 
observations by the XMM-Newton observatory (lansen et al. 2001), 
that provide the most detailed view of the X-ray spectrum of this 
ULX, and hence the best constraints upon the accretor, to date. 



2 THE DATA 

The new XMM-Newton data were obtained in two observations sep- 
arated by a gap of 109 days. We present details of the observa- 
tions in Table Q All three EPIC detectors were operated in full 
frame imaging mode with the medium filter in place. Data for the 
analysis were extracted from the pipeline product event lists us- 
ing the appropriate XMM-Newton SAS v.5.4.1 tools. The pn data 
were filtered to leave only events with flag = and pattern 4; 
MOS data were filtered for pattern ^ 12 and the "Jxmmea.em" 
flag. We initially created images from these events in the 0.3 - 10 
keV range, and checked the background for flaring by producing 
a light curve from the whole detector excluding the target source. 
The background count rate in the first observation was consistently 
low (< 10 count in the pn detector), so no time filters were ap- 
plied. However, the second observation was heavily contaminated 
by background flaring, so time intervals during which the back- 
ground rate in the pn exceeded 30 count s^^ (in 100 s bins) were 
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Figure 1. (top) The central regions of the combined XMM-Newton EPIC 
image of NGC 5204. The image is the sum of the data from all three detec- 
tors and both epochs (subject to the event selections discussed in the text). 
An approximate exposure correction has been implemented. The image size 
is an 8 X 8 arcminute^ box, the data are convolved with a 1 pixel (= 4 arc- 
sec) HWHM 2-D Gaussian mask, and the intensity logarithmically-scaled, 
for illustrative purposes. Contours are overlayed at 5 and 10 count pixel"^ . 
The data extraction region for NGC 5204 X-1 is shown by the solid cir- 
cle; background extractions were made in the dashed (first observation) and 
dot-dashed (second observation) circles respectively, (bottom) The corre- 
sponding OM UVWl data from the first observation. The data have been 
astrometrically corrected using reference stars. The position of NGC 5204 
X-1 is highlighted by the 10-arcsec radius circle. 



excluded. This resulted in the loss of ~ 50% of the data. The re- 
maining exposure times (as listed in the GTI file extensions for the 
on-axis CCDs, which are also corrected for telemetry drop-outs) 
are given in Tabled The filtered and co-added 0.3 - 10 keV im- 
ages are displayed in Figure Q Contemporaneous OM data (in the 
UVWl, UVM2 and UVW2 filters) were also obtained. The UVWl 
image of NGC 5204 is displayed in FigureQfor comparison with 
the EPIC image (although we defer a full analysis of the OM data 
to a later work). 

NGC 5204 X-1 was detected at its expected position 
(13''29'"38.6", +58°25'06" (J2000); see e.g. Roberts et al. 2004). 
Source spectra and lightcurves were extracted in a 9-pixel (= 36 
arcsec) radius circle centred on the source, chosen to sit inside the 
chip gaps close to the pn aimpoint. An 11 -pixel radius region cen- 



Table 1. Details oi the XMM-Newton observations. 



Observation ID 


Date 


Useful exposure (s) 




(yyyy-mm-dd) 


pn 


MOS 


0142770101 


2003-01-06 


17047 


18665 


0142770301 


2003-04-25 


8609 


8645 



tred ~ 2 arcmin from NGC 5204 X-1, on the same pn chip and at 
approximately the same distance from the readout nodes, was used 
to produce background data (see Figure0. The spectral products, 
including RMF and ARF matrices (which contain the appropriate 
corrections for the source extraction region size), were created us- 
ing the ESPECGET tool, and the source spectra were binned to a 
minimum of 20 counts per bin before analysis in XSPEC vll.3. 



3 RESULTS 

Before embarking upon a detailed analysis of the source properties, 
it is worth clarifying whether the source signal is contaminated by 
other X-ray emitters within NGC 5204. In particular, Chandra ob- 
servations of NGC 5204 X-1 highlight a second point-like source 
located ~ 15 arcsec to the south-east of NGC 5204 X-1 (Roberts 
et al. 2001; Roberts et al. 2004), well within the data extraction 
region used in this analysis. If this source maintains its historical 
flux level with respect to NGC 5204 X-1, where its Chandra count 
rate was typically <C 2% of that of the ULX, it will not signif- 
icantly contaminate the data. We checked whether this is true by 
extracting the radial profile of NGC 5204 X- 1 from the combined 
EPIC image. If the fainter source had brightened considerably then, 
since it is located ~ 15 arcsec from X-1, it should be marginally 
resolved in the data and would cause an excess above the expected 
XMM-Newton point spread function (PSF) in one or more radial 
bins. However, the off-peak drop-off in the PSF was smooth, im- 
plying that contamination from the fainter source is not a strong 
effect. We also note that the apparently extended emission outside 
the source extraction region (c.f. Figure can adequately be ex- 
plained by the ~ 15% of the source flux in the wings of the PSF 
that falls outside the extraction region^ . Hence there is no evidence 
for any strong contamination of the NGC 5204 X-1 signal from 
other X-ray sources in NGC 5204. 

We tested for temporal variability from NGC 5204 X-1 by ex- 
tracting light curves from both sets of data. The mean 0.3 - 10 
keV count rate from the combined EPIC detectors rose by ~ 34% 
between the two observations (from 0.86 to 1.15 count s"'^), ex- 
tending the historical pattern of low-amplitude long-term variabil- 
ity from this ULX (c.f. Roberts et al. 2004, Figure 7). Using the 
spectral models discussed below gives 0.5 - 2 keV fluxes of ^ 0.74 
and 1.05 x 10~^^ erg cm~^ s~^, which sit within the previously- 
observed flux range. However, neither observation showed strong 
evidence for short-term variability. For example, extracting a power 
density spectrum in the full 0.3 - 10 keV band from the longer first 
observation reveals only a slight excess of power above the pois- 
son noise level at low frequencies (< 0.01 Hz). Binning the data 
to 100-s intervals allowed this variability to be detected as excess 

^ The apparent north-south structure in FigurelTI mirroring the morphology 
of the galaxy, may in part be due to the presence of pn chip gaps. However, 
some of the extension on the periphery of the contour may be due to real 
faint X-ray emitters in NGC 5204. 
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Table 2. Spectral fits. 



Model parameter 



2003 January 6 



2003 April 25 



Model 1: WA 

Nh'' 

r 

Model 2: WA 

Nh 



log (/) 9 

^BMC " 

Model 3: WA 

Nh 

r 

xVdof 
Model 4: WA 

r 

Ad 



PO ' 



0.56 ± 0.08 



2.14 



+0.05 



2+0.14 
-0.13 



-0.04 

(4.15 ±0.16) X 10- 
573.6/540 

BMC 

0/ 

0.16 ± 0.01 
0.96 ± 0.06 
0.48^ 
(1.1 ±0.1) X 10"^ 
530.6/539 
(DISKBB + PO) 

0.21 ± 0.03 

1 7+36 
-^'-10 

1.97 ±0.07 

+°:^) X 10-* 

530.9/538 



1 Ol+oi"^ 

^■^^-0.12 
^■■^^-0.06 

(6.76«:^3^) X 10- 
444.6/394 



Q4j^+0.18 

0.17 ±0.02 



1 oq+°-^2 
^■'J^-o.ii 

,+0.15 
-0.14 



0.16"* 



(1.8toj) X 10- 
380.6/392 



(3.4 



:,+0.34 
-0.29 
1+0.04 
-0.03 
^+88 
-28 
□+0.12 
-0.13 



1.22: 
0.21^ 
45^ 
2.09"* 
(4.7+0-;) X 10- 
387.3/392 



(PO 



+ DISKBB) 
+0.48 
0.28 
+0.43 
0.26 

(4.7l?,:^) X 



1.68: 

3.34"* 



-0.4' 
2.15 



10- 



n+0.45 
-0.34 
^+0.41 

-0.33 
^+10) X 10-4 



2.00: 

3.34"* 



+0.23 
-0.25 

{2.2tll) X 10- 
526.5/538 

1.6 X 10-12 
4.4 X 10^9 



2.78 



+ 1.65 
-0.66 



(8.41^^1'') X 10- 
409.6/392 



2.0 X 10-12 
5.5 X 10^9 



Notes: " Models are abbreviated to XSPEC syntax. The components are: 
WA - cold absorption; PO - power-law continuum; BMC - bulk motion 
Comptonization spectrum (see text); DISKBB - MCDBB. Absorption 
column external to our Galaxy in units of lO^^ atom cm~^. A Galactic 
absorption column of 1.5 X 10^*^ atom cm— ^ is assumed (see text). 
XSPEC model normalisations. Disc colour temperature in units of keV. 
" Power-law spectral index, f Value fixed at that shown. ^ Illumination 
parameter (see e.g. Borozdin et al. 1999). ^ Inner disc temperature in keV. 
* Flux from Model 4 (which agrees with Models 2 and 3 to ~ 3%), in 
units of ergcm-2 g- ^ (0.3 - 10 keV). ^ Derived X-ray luminosity, in 
erg s~i, for a distance of 4.8 Mpc. 



variance, albeit at a very low amplitude (~ 3% of the total signal). 
The second observation showed no similar variability, with a 95% 
upper-limit to its excess variance (on the same timescale) of « 2% 
of the total signal. The following analysis therefore focuses solely 
on the time-averaged spectrum extracted from each observation. 

The spectra from all three EPIC detectors were fit concur- 
rently in the 0.3 - 10 keV range. We allowed for differences be- 
tween the detectors by including a constant multiplicative factor 
in each model, that was fixed at unity for the pn data and allowed 
to vary for the two MOS detectors. Typical values obtained were 
0.88/0.86 and 0.93/0.88 in the first/second observations for MOSl 
and MOS2 respectively, implying a disagreement of < 15% in the 
absolute flux measurements from the instruments (we quote an av- 
erage of the three measurements where we give fluxes). We applied 
a fixed column of 1.5 x 10^'' atom cm-^ to each spectral model to 
represent the foreground column through our galaxy (interpolated 
from Stark et al. 1992), plus a second absorption component left 
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Figure 2. Residuals to the power-law continuum fits in both epochs. The 
spectral data are rebinned to aid clarity, and data from all three EPIC de- 
tectors are shown with MOS1/MOS2 data points identified by open cir- 
cles/squares. 



free to fit the data, which represents additional absorption within 
NGC 5204 and/or intrinsic to the ULX. 

The spectra appeared relatively smooth and featureless, hence 
we began by fitting a number of simple continuum models to the 
data. First we tried a power-law continuum model; this gave rea- 
sonable fits, with xt ~ 1-06 and ~ 1.12 in the first and second 
observations respectively. The details of these fits are given in Ta- 
ble|2|(Model 1). We note that the overall spectrum of NGC 5204 X- 
1 appears to have hardened substantially in the two years between 
the Chandra and XMM-Newton observations, from F ~ 2.8 — 3 to 
~ 2.15 — 2.35 (cf. Roberts et al. 2004). However, Figure|2|demon- 
strates that the same substantial systematic residuals to the power- 
law fit were present in both XMM-Newton observations, implying 
that more complex models are required to fit the data. In partic- 
ular the spectral shape is inconsistent with any single component 
model with a convex continuum form, such as a classical black- 
body, a thermal bremsstrahlung model, or the "multi-colour disc 
blackbody" spectrum used to represent the optically-thick thermal 
emission from an accretion disc (hereafter MCDBB; see e.g. Mit- 
suda et al. 1984). As expected the spectral fits of these "single com- 
ponent" models were rejected at a high statistical significance in 
both sets of data (xl > 2 for the blackbody and MCDBB, and 
xi ^ 1-4 for the thermal bremsstrahlung model). 

We next considered models previously employed to explain 
the spectra of ULXs. A broken power-law continuum fit (e.g. 
Kaaret et al. 2003) could represent X-ray emission produced in the 
Compton-scattering of photons from the accretion disc and/or the 
companion star by a relativistic jet (cf. Georganopoulos, Aharo- 
nian & Kirk 2002). However, we ruled this model out as it predicts a 
spectral downturn at low energies (i.e. TiowE < ^MghE), whereas 
both our spectra show an upturn at '--^ 1 keV (cf. Figure|2j. A sec- 
ond possible model was a Comptonized disc model. Kubota, Done 
& Makishima (2003) suggest that ULXs with a power-law-like X- 
ray spectrum are actually in a state analogous to the very high state 
of Galactic black hole binaries, where their emission is dominated 
by a strongly-Comptonized accretion disc. We first approximated 
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this spectrum using tlie "COMPTT" model of Titarcliuk; (1994) in 
XSPEC, and found better fits than the power-law continuum model 
in both cases {xl ~ 1). However, the absence of data above 10 
keV leaves the parameters of this model poorly constrained - rea- 
sonable values (including a cool 130 - 140 eV seed photon temper- 
ature) were only obtained after fixing the electron temperature of 
the scattering medium at 50 keV. As a second approximation of a 
Comptonized disc spectrum we used the bulk motion Comptoniza- 
tion model ("BMC", see Schrader & Titarchuk 1999 and references 
therein), which has the twin advantages that it does not require hard 
(> 10 keV) X-ray data to constrain its fitting, and it provides an es- 
timate of the inner accretion disc temperature. This model again 
provides a very good fit to each dataset (Model 2 in Table|2}. Most 
importantly, it suggests that a "cool" accretion disc, with a temper- 
ature of ~ 160 - 170 eV, underlies the spectrum of NGC 5204 X- 1 . 
This is consistent with the presence of an accreting IMBH. 

To confirm this result, and place it in the context of other 
recent ULX observations, we then fit the data with a combined 
MCDBB -I- power-law model. This empirically describes the spec- 
tra of black hole X-ray binaries (BHBs; cf. McClintock & Remil- 
lard 2004; hereafter MR04), and has been used to demonstrate ev- 
idence for cool accretion discs in ULXs (e.g. Miller et al. 2003)^. 
We find, as expected given the Comptonization fits, that the data are 
very well fit by this "standard" model, displaying a cool (fcT ~ 0.2 
keV) MCDBB component plus a hard power-law continuum. The 
best-fit parameters are listed in Table|2|as Model 3. In each observa- 
tion the detection of the "cool" disc component is apparently very 
robust, with F-test significances of > 6cr and > 7a respectively 
when compared to the simple power-law continuum fit. 

Stobbart et al. (2004) have recently demonstrated that a vari- 
ant of the standard model provides the best empirical descrip- 
tion of the spectrum and spectral changes exhibited by a marginal 
(Lx ^ 10^® erg s"^) ULX in the very nearby galaxy NGC 55. In 
their analysis, the ULX spectrum changes over the course of two 
contiguous 30-ks observations, initially having a soft power-law 
continuum form (F ~ 4). As the observation progressed a sec- 
ond, curved continuum spectral component was required to model 
an increase in the flux predominantly above 1 keV, that was well- 
fitted by a kTi„ ~ 0.9 keV MCDBB spectrum. Whilst utilis- 
ing the same pair of components as in the standard description, 
this particular model is distinguished by the power-law continuum 
dominating at low energies, the reverse of the usual prescription. 
[A caveat, of course, is that this strictly applies only to observa- 
tions carried out in the limited 0.3 - 10 keV bandpass of XMM- 
Newton.] This "non-standard" spectral model was also indepen- 
dently established to provide a good fit for the nearest persistent 
extra-galactic ULX, M33 X-8, with F ~ 2.5 and kTi,^ ~ 1.16 keV, 
by Foschini et al. (2004). The NGC 55 and M33 X-8 ULXs have 
Lx ~ 1 — 2 X 10"^^ erg and, therefore, are both at the low- 
luminosity end of the ULX regime. The question of whether the 
non-standard model provides at least an empirical representation 
of the spectra of more luminous ULXs therefore merits some con- 
sideration. On this basis we have attempted to fit the non-standard 
model to the current data. In fact an acceptable fit to the data from 
both epochs was obtained by this approach (Model 4 in Table|2j. 
We explicitly demonstrate that the data allows two minima - cor- 



^ This description is in essence a simplified version of the BMC model, 
in that it employs the same spectral components, but it does not treat their 
relative fluxes in a physically self-consistent manner, or truncate the power- 
law continuum at low energies. 
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Figure 3. The minima in x^-space, con'esponding to Models 3 & 4, shown 
as a function of fcTi,j and F. The contours are displayed at Ax^ = 27, 9 
and 3 above the best fit for the observation (with the largest deviation 
mapped by the outermost contour). 



responding to Models 3 & 4 - by plotting the minima in kTin — F 
confidence space in Figure|3| 

In Figure |4| we show the data and best-fit model (Model 4) 
for the first observation (the fit quality in the second observation 
is consistent with that shown), and include a direct comparison of 
the underlying forms of Models 3 & 4. The composite spectrum is 
almost identical in the two cases, with the largest divergence evi- 
dent at high energies (> 3 keV), where Model 4 shows a spectral 
downturn due to the distinct curvature in the dominant MCDBB 
component. We have further investigated whether the data best sup- 
ports Model 3 or Model 4 by explictly searching for evidence of 
this high-energy break in the spectral data. We restricted our anal- 
ysis to the data points above 2 keV, and modelled the spectra using 
both single and broken power-law continua. We find that a broken 
power-law is not required by the second (April 25) dataset, which 
is very well-modelled above 2 keV by a F ^ 2.1 power-law contin- 
uum (x^/dof = 84.9/83), consistent with Model 3. However, whilst 
the first (January 6) dataset is also well-fitted by a F ~ 2.45 power- 
law continuum above 2 keV (x^/dof = 137.4/136), we find that the 
fit is improved at the ~ 3cr confidence level (according to the F- 
statistic, i.e. Ax^ = 9.7 for two additional degrees of freedom) by 
a broken power-law fit with Fi = 1. 8^0.3, ^break = 4.91q;4 keV, 
and F2 = S.l^Q g. Hence the first observation shows marginal ev- 
idence for the presence of a spectral downturn at 5 keV, which (in 
combination with the slightly better fit to the full spectrum) favours 
Model 4 as the best description of its spectral shape. 
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2003 January 6, data and Model 4 
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Figure 4. Best fit spectra to the first observation using the variants of 
the MCDBB + power-law continuum model, (top panel) The data (with 
MOS1/MOS2 points distinguished by open circles/squares as before) and 
best fit Model 4 spectrum (solid line) are shown in the upper part of the 

panel, and the fit residuals in Ax space in the lower part, (centre & bottom 
panels) A comparison of the variants of the MCDBB + power-law contin- 
uum models. The overall fit is shown as a solid line, with the contributions 
from the MCDBB and power-law continuum distinguished as dotted and 
dashed Unes respectively. 



4 DISCUSSION 
4.1 IMBH models 

Both the Comptonization models and the accepted empirical model 
for BHB spectra point towards the presence of a cool accretion disc 
signature in the X-ray spectrum of NGC 5204 X-1. As discussed 
by previous authors (e.g. Miller et al. 2003), this suggests that an 
IMBH is present in this system. Indeed, using equations (5) and 
(8) of Makishima et al. (2000), the best-fit MCDBB parameters 
translate into masses of 256(cos(i))~ 5 and 412(cos(i))~2 M© 
(where i is the inclination of the accretion disc with respect to the 



observer) for the first and second observations respectively. Clearly 
this places the accreting object well into the IMBH regime^. 

The optical counterpart of NGC 5204 X-1 (HST-3 of Goad et 
al. 2002) lies in close proximity to a probable young stellar cluster, 
HST-1. Liu et al. (2004) resolve HST-1 into a point-like red compo- 
nent (a possible yotmg cluster of a few thousand stars) and a more 
diffuse blue component, probably an OB-association. This prox- 
imity is potentially important given recent results suggesting that 
IMBHs can be formed by runaway stellar collisions in young, dense 
stellar clusters (Portegies-Zwart et al. 2004 and references therein). 
The projected separation of HST-3 from HST-1 is only ~ 26 pc; if 
the ULX originated in the cluster, this distance must have been tra- 
versed in less than the ^ 10 year lifetime of the BO lb secondary. 
This is possible with a transverse velocity as small as ~ 2 km 
for the ULX relative to the stellar cluster. If we consider that stellar- 
mass black holes in our Galaxy are seen with spatial velocities, pre- 
sumably imparted by a kick at the time of the black hole formation, 
of the order ~ 100 km s^^, then the minimum velocity require- 
ment appears entirely plausible for an ULX harbouring a ~ 500 
Mq IMBH, formed in, and subsequently kicked-out of the clus- 
ter (cf. Zezas et al. 2002 for similar calculations for ULXs in the 
Antennae). 

However, there is an anomaly in the parameterization of spec- 
tral Model 3, from which the IMBH result is derived. The detection 
of an IMBH is based on a reasonable extrapolation of stellar-mass 
black hole properties, as represented by the empirical MCDBB + 
power-law spectral model, to the IMBH regime, with the "cool" 
accretion disc signature the key evidence for an IMBH. Because 
of the extreme luminosities of ULXs it is generally assumed that 
the accretion disc is in a "high" state, where it extends to the in- 
nermost circular orbit around the black hole. This must correspond 
to the classic "high/soft" state in a stellar-mass BHB (also known 
as the "thermal dominated" state, cf. MR04). However, in BHBs 
the MCDBB in this state is typically accompanied by a power- 
law component with a measured photon index > 2.4. This is also 
true for the "very high" (steep power-law) state. In fact, such steep 
slopes are theoretically expected from Comptonized accretion disc 
spectra of sources close to the Eddington limit (Chakrabarti & 
Titarchuk 1995). Assuming that the XMM-Newton passband allows 
a reasonable measurement of the power-law slope, then such steep 
slopes are not present in our data when parameterized by Model 3. 
This anomaly is actually seen in several ULXs with inferred soft ac- 
cretion discs (e.g. Miller et al. 2003; Miller, Fabian & Miller 2004), 
though not all of them (cf. Soria et al. 2004; Dewangen et al. 2004). 
The anomalous power-law slopes derived for ULXs with the IMBH 
model are generally < 2, which is actually more consistent with the 
"low/hard" state in BHBs. 

So could some ULXs instead be low/hard state sources? This 
state generally shows much cooler accretion disc temperatures than 
the high/soft state in BHBs (possibly due to the truncation of the 
inner edge of the accretion disc at large radii from the black hole; 
MR04), with observed values of kTin similar to the temperatures 
derived in the IMBH model. Hence a low temperature disc is not 
necessarily an unambiguous indicator of the presence of an IMBH. 
[We note that an IMBH in the low/hard state would, by extended 
analogy, have an even cooler disc than that observed {kTi„ ~ 10s 
of eV).] However, the luminosities associated with the low/hard 



3 In this calculation we assume a = 1 in equation (8) of Makishima et al. 
(2000), i.e. the Schwartzschild metric. A Kerr metric would imply a < 1 
hence a larger black hole mass. 
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state in BHBs are generally very low by comparison to ULXs 
(Lx ~ 10^® erg s~^). This characteristic lack of X-ray luminos- 
ity appears to rule out stellar-mass black holes in the low/hard state 
as candidates to explain the power-law slope of Model 3. It seems 
therefore that the classic BHB states may not always be exact analo- 
gies of ULX behaviour, which could be problematic for the IMBH 
interpretation. 

4.2 "Non-standard" models 

An alternative description of the X-ray spectrum of NGC 5204 X-1 
is provided by Model 4, which as previously stated has been shown 
to provide the best-fit model to high-quality XMM-Newton spectra 
in two ULXs (Stobbart et al. 2004; Foschini et al. 2004). In this 
description, both the inner-disc temperature and power-law slope 
are well within the range of known behaviour for the high/soft and 
very high states in Galactic stellar-mass BHBs (fcTj„ > 0.7 keV; 
r > 2.5)". Whilst the inferred disc temperatures for NGC 5204 
X-1 in this model do appear high for a stellar-mass system, sim- 
ilar temperatures have been observed from BHBs in the high/soft 
and very high states (MR04) and, interestingly, a high disc temper- 
ature is consistent with the "slim disc" regime for near- and slightly 
super-Eddington accretion (Watarai, Mizuno & Minishige 2001; 
Ebisawa et al. 2003). A stellar-mass system is also not ruled out 
from the proximity of the ULX to a possible young stellar cluster. 
Indeed, a very plausible scenario is that the ULX is simply a high- 
mass X-ray binary recently formed within, and kicked-out of, the 
OB-association. 

However, there are physical problems with the non-standard 
model. The MCDBB normalisation translates into inner disc radii 
(Rin) of ~ 26.5 and ~ 16 km. These are very problematic as 
they are somewhat smaller than the 30 km Schwarzschild radius of 
a IO-M0 black hole. However, comparable discrepancies are also 
seen in some Galactic systems. For example, Ueda et al. (2002) 
use ASCA data to measure Ri„ ~ 30 km for two observations of 
GRS 1915+105, a black hole that has a dynamically-determined 
mass of 14 ± 4 M© (Greiner, Cuby & McCaughrean 2001). Again, 
apparently small inner-disc radii are also a prediction of the slim 
disc model (Watarai, Mizuno & Minishige 2001). We also note the 
apparent small inner-disc radius in NGC 5204 X-1 could in part 
be explained by an inclination i > which will serve to increase 
Rin- The masses derived from the MCDBB model normalization 
are within the stellar- mass black hole regime, at 3(cos(i))~2 and 
1.8(cos(i))~2 Mq for the two observations. This suggests that 
NGC 5204 X-1 could contain a stellar-mass black hole. For a ~ 10 
M0 black hole, this would require a radiation rate at ~ 4 — 5 
times the Eddington luminosity. We note this is well within the 
predictions of Begelman (2002) and King et al. (2001) for truly 
super-Eddington/apparently super-Eddington accretion discs, and 
such luminosities have been observed (albeit briefly) from Galactic 
stellar-mass systems (MR04). 

The most severe problem with the non-standard spectral de- 
scription is the domination of the power-law continuum at low en- 
ergies. This spectral form has not (to the best of our knowledge) 
previously been observed for any Galactic BHB. Furthermore, it 
is physically inconsistent with the standard accretion disc plus 
Comptonized-corona model for BHBs, as one would not expect to 

* Indeed, the derived values in observation 1 are extraordinarily close to 
the fcTi„ ~ 2.2 keV, T ~ 3.5 fit quoted for GRS 1915+105 in its high/soft 
state by MR04. 



see the power-law continuum extend below the peak emissivity of 
the MCDBB component, due to a lack of seed photons available to 
be Compton-scattered. It is very unlikely that the corona is simply 
seeing an alternate source of seed photons; for instance it would 
require a radically different disc emission profile for these photons 
to originate in the cool outer regions of the accretion disc. It is also 
very unlikely that the secondary star could seed this component; as- 
suming equal primary and secondary masses of ~ 10 Mq, with the 
BO lb secondary overflowing its Roche lobe (cf. Liu et al. 2004), 
and a corona extending out to 100 Schwarzschild radii, then the 
coronal region intercepts less than one in a million photons emit- 
ted by the star. As the power-law component emits ~ lO**^ photon 
s~^ above 0.3 keV (after correcting for absorption) this requires 
implausibly high stellar fluxes. 

So can the dominant soft power-law continuum instead be ex- 
plained by an alternative physical process? Foschini et al. (2004) 
discuss two possibilities for M33 X-8; mildly-relativistic jets, and 
an outflowing wind. A mildly-relativistic jet is certainly a possibil- 
ity in M33 X-8 due to the detection of a radio counterpart. How- 
ever, Galactic black holes tend to display spectrally-hard power-law 
continua (F ~ 1.5 — 2) during episodes of radio jet emission, pre- 
dominantly when in the low/hard spectral state (though jets can 
also be seen from very high state sources, e.g. GRS 1915+105; 
Reig, Belloni & van der Klis 2003). Such a hard power-law con- 
tinuum is inconsistent with the X-ray spectral forms seen in these 
ULXs. The alternative, an outflowing wind, may be likely if ULXs 
are powered by accretion from a high-mass star onto a stellar-mass 
black hole. For instance, Podsiadlowski, Rappaport & Han (2003) 
demonstrate that the rate of matter transfer from the star could po- 
tentially be far higher than that required to reach the Eddington 
limit in these systems. Presumably this excess material must either 
be accreted onto the black hole in a radiatively-inefficient manner, 
or more likely discharged from the system in an outflow (see Begel- 
man 2002, King 2002). This scenario may be particularly appropri- 
ate for NGC 5204 X-1, which possesses a BO lb star companion 
(Liu et al. 2004). However, whilst such a wind could emit a soft 
X-ray spectrum (e.g. King & Pounds 2003), it should have a ther- 
mal blackbody spectrum. Hence this also appears to be unable to 
provide an origin for the dominant soft power-law. 

The lack of an obvious explanation for the origin of the soft 
power-law continuum seems to provide insurmountable physical 
problems for this model. However, given that it provides a clear 
description of the shape of the X-ray spectra of M33 X-8 and the 
NGC 55 ULX (and is marginally preferred in the first observation 
of NGC 5204 X-1), Model 4 must demonstrate at least subtle differ- 
ences to, or perhaps greater complexity than, the physical processes 
underlying in the standard accretion disc plus optically-thin corona 
model. Wc speculate that one physical picture that may explain this 
distinct spectral shape is presented by Zhang et al. (2000). They 
describe accretion discs around stellar-mass black holes by a three- 
layer model roughly analogous to the Solar corona, composed of 
a cold (0.2 - 0.5 keV) inner-disc layer, a warm and optically-thick 
atmosphere to the disc (r ~ 10, kT ~ 1 — 1.5 keV) and a much 
hotter optically-thin corona. We note that such a model will in- 
evitably produce a Comptonized component down to the energies 
suggested by the soft power-law continuum in Model 4. Crucially 
the warm, optically-thick accretion disc atmosphere will also pro- 
duce a spectral turnover at energies of a few keV, similar to that we 
have modelled by a MCDBB. It is therefore possible that at least 
some ULXs could have X-ray spectra dominated by this optically- 
thick warm accretion disc atmosphere. This will be investigated in 
detail in future work. 
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4.3 Concluding comments 

These observations imply a degree of spectral ambiguity as mod- 
els supporting the presence of an IMBH (Models 2 and 3 of Ta- 
ble|2j and the non-standard alternative (Model 4) provide equally 
acceptable fits to the available data. In our view the balance of the 
argument still appears to favour the presence of an IMBH in NGC 
5204 X-1, although its hard power-law slope remains as somewhat 
of an anomaly. However the contrasting interpretation of the data 
discussed in this paper clearly merits further investigation. 
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